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I. INTRODUCTION 



Ying Fan0 Yan-Qing MaQ and Kuang-Ta ChacH 
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Relativistic corrections to the color-singlet J/tp hadroproduction at the Tevatron and LHC are 
calculated up to 0{v'^) in nonrelativistic QCD (NRQCD). The short distance coefficients are ob- 
tained by matching full QCD with NRQCD results for the subprocess g + g ^ J/ip + g- The long 
distance matrix elements are extracted from observed J/ip hadronic and leptonic decay widths up 
to 0{v^). Using the CTEQ6 parton distribution functions, we calculate the LO production cross 
. . . sections and relativistic corrections for the process p + p{p) J/tp + X at the Tevatron and LHC. 

On ' We find that the enhancement of 0{v^) relativistic corrections to the cross sections over a wide 

, range of large transverse momentum pt is negligible, only at a level of about 1%. This tiny effect is 

^ ' due to the smallness of the correction to short distance coefficients and the suppression from long 

distance matrix elements. These results indicate that relativistic corrections can not help to resolve 
$-H , the large discrepancy between leading order prediction and experimental data for J/tp production 

at the Tevatron. 

< 

(N 

^! Nonrelativistic QCD (NRQCD) [H is an effective field theory to describe production and decay of heavy quarkonium. 
(-H ' In this formalism, inclusive production cross sections and decay widths of charmonium and bottomnium can be factored 
into short distance coefficients, indicating the creation or annihilation of a heavy quark pair, and long distance matrix 
elements, representing the evolvement of a free quark pair into a bound state. The short distance part can be calculated 
perturbatively in powers of coupling constant as, while the nonpcrturbative matrix elements, which are scaled as v, 
the typical velocity of heavy quark or antiquark in the meson rest frame, can be estimated by nonpcrturbative methods 
or models, or extracted from experimental data. 
I One important aspect of NRQCD is the introduction of color-octet mechanism, which allows the intermediate 
■rj" ' heavy quark pair to exist in a color-octet state at short distances and evolve into the color-singlet bound state at 
t;-^ . long distances. This mechanism has been applied successfully to absorb the infrared divergences in P-wave[ll [3, Q 
and D-wave0,Q decay widths of heavy quarkonia. In Ref. the color-octet mechanism was introduced to account 
On . for the J/tP production at the Tevatron, and the theoretical prediction of production rate fits well with experimental 
' data. However, the color-octet gluon fragmentation predicts that the JJip is transversely polarized at large transverse 
^ momenturri pt, which is in contradiction with the experimental data[7|. (For a review of these issues, one can refer 
"'H to Refs.jl, [§, [lO'l). Moreover, in Rcfs.fTll. [T^ it was pointed out that the color-octet long distance matrix elements 
K% of J/ip production may be much smaller than previously expected, and accordingly this may reduce the color-octet 
contributions to J/ip production at the Tevatron. 

In the past a couple of years, in order to resolve the large discrepancy between the color-singlet leading-order (LO) 
predictions and experimental measurements [isi fl3. [isj of J/ip production at the Tevatron, the next-to-leading-order 
(NLO) QCD corrections to this process have been performed, and a large enhancement of an order of magnitude 
for the cross section at large pt is found[l^, But this still can not make up the large discrepancy between the 
color-singlet contribution and data. Similarly, the observed double charmonium production cross sections in e~^e~ 
annihilation at B factories [Tsl IToj also significantly differ from LO theoretical predictions [1^ . Much work has been 
done and it seems that those discrepancies could be resolved by including NLO QCD corrections [2ll [2^. [23l. [2^ and 
relativistic corrections [l^ |2^. One may wonder if the relativistic correction could also play a role to some extent in 
resolving the long standing puzzle of J/tp production at the Tevatron. 

In this paper we will estimate the effect of relativistic corrections to the color-singlet J/ip production based on 
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NRQCD. The relativistic effects are characterized by the relative velocity v with which the heavy quark or antiquark 
moves in the quarkonium rest frame. According to the velocity scaling rules of NRQCD [27j. the matrix elements 
of operators can be organized into a hierarchy in the small parameter v. We calculate the short distance part 
perturbatively up to 0{v'^). In order to avoid model dependence in determining the long distance matrix elements, we 
extract the matrix elements of up to dimension-8 four fermion operators from observed decay rates of J/'^[28|. We find 
that the relativistic effect on the color-singlet J/tp production at both the Tevatron and LHC is tiny and negligible, 
and relativistic corrections can not offer much help to resolve the puzzle associated with charmonium production at 
the Tevatron, and other mechanisms should be investigated to clarify the problem. 

The rest of the paper is organized as follows. In Section II, the NRQCD factorization formalism and matching 
condition of full QCD and NRQCD effective field theory at long distances are described briefly, and then detailed 
calculations are given, including the perturbative calculation of short-distance coefficient, the long-distance matrix 
elements extracted from experimental data, and the parton level differential cross section convolution with the parton 
distribution functions (PDF). In Section III, numerical results of differential cross sections over transverse momentum 
pt at the Tevatron and LHC are given and discussions are made for the enhancement effects of relativistic corrections. 
Finally the summary of this work is presented. 



II. PRODUCTION CROSS SECTION IN NRQCD FACTORIZATION 

According to NRQCD factorization[l], the inclusive cross section for the hadroproduction of J/tp can be written as 

^(g + 5 - J/V^ + .9) = E ^(0|O,{/'/l0). (1) 

n ^ 

The short-distance coefficients Fn describe the production of a heavy quark pair QQ from the gluons, which come from 
the initial state hadrons, and usually expressed in kinematic invariants, rric is the mass of charm quark. The long- 
distance matrix elements (0|C',i^'^|0) with mass dimension dn describe the evolution of QQ into J/ijj. The subscript 
n represents the configuration in which the cc pair can be for the J/ip Fock state expansion, and it is usually denoted 
as 71 = '^^^^L^j'^\ Here, S,L and J stand for spin, orbital and total angular momentum of the heavy quarkonimn 
respectively. Superscript 1 or 8 means the color-singlet or color-octet state. 

For the color-singlet '^S'l cc production, there are only two matrix elements contributing up to 0{v'^): the leading 

order term {0\O''^^ {'^ s\^^)\0) and the relativistic correction term {0\V''^^(^S^i^)\0). Therefore the differential cross 
section takes the following form 



^(g + g^J/^ + g) = ^l!^ (0|O'V^ ^5 ^ ) |0) + (0|7'^/'^f ) |0) + 



(2) 



And the explicit expressions of the matrix elements are[lj: 

^O|0.//V(3^[i])|O^ = (0|xVV(a;av,)VV^x|0), 

{0\V'/^eS^^^m = (0|l[xVV(a;a^)7^V^(-^D)2x + /i.c.]|0), (3) 

where ip annihilates a heavy quark, x creates a heavy antiquark, and are operators creating and annihilating 
J/tp in the final state, and D = D — D 

In order to dete 
NRQCD is needed: 



In order to determine the short distance coefficients and G('^S'[^'), matching condition of full QCD and 



dt 



pert QCD 



pert NRQCD 



The differential cross section for the production of charmonium J/ip on the left hand side of Eq. (U can be calculated 
in perturbative QCD. On the right hand side the long distance matrix elements are extracted from experimental data. 
Quantities on both sides of the equation are expanded at leading order of as and next-to-leading order of w^. Then 
the short distance coefficients F(^S'P') and G'('^S'P') can be obtained by comparing the terms with powers of on 
both sides. 
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FIG. 1: Typical Feynman diagram for ''sf ' cc hadroproduction at LO. 



A. Perturbative Short Distance Coefficients 



We now present the calculation of relativistic correction to the process g + g ^ J/4> + g- In order to determine the 
0{v^) contribution in Eq.®, the differential cross section on the left hand side of Eq.(|3]) or equivalently the QCD 
amplitude should be expanded up to 0{v'^). We use FeynArts [s^l to generate Feynman diagrams and amphtudes, 
FeynCalc [si'] to handle amplitudes, and Fortran to evaluate the phase space integrations. A typical Feynman diagram 
for the process is shown in Fig. [TJ The momenta of quark and antiquark in the lab frame are(26» .29.. ,32j: 

ip + g = L{^Pr + qr), 

\P-1 = L{\Pr-qr), (5) 

where L is the Lorenz boost matrix from the rest frame of the J/V' to the frame in which it moves with four momentum 
P. Pr — {'2Eq, 0), Eq — VTO^ + I^, and 2qr = 2(0, q) is the relative momentum between heavy quark and antiquark 
in the J /ip rest frame. The differential cross section on the left hand side of Eq. ^ is 



da , 



{g + g^J/^j + g) ^ ^y^^M[g + g ^ J/^ + g)\\Q\0''^{'S^^^m, (6) 



pert QCD IQ'Ks'^ ^ 



where (0|C''^/'^('^S'[^')|0) is the matrix element evaluated at tree level, and summation/ average of the color and po- 
larization degrees of freedom for the final/initial state has been implied by the symbol The amplitude for the 
color-singlet process g{p\) + g{p2) J/ip{P3 — P) + 9{P4) is 



/ZTr[C[iln(i)X 



M{g + g^ J/^j + g) = . —Trid^W^^M^""], (7) 



where denotes the parton-level amplitude amputated of the heavy quark spinors. The factor y "g- comes from 

the normalization of the composite state pS'|^')i5i]. Here the covariant projection operator methodjH, [131 is adopted. 
For a color-singlet state, the color projector C'^' = -^=- The covariant spin-triplet projector 11'^^) in ([7]) is defined by 

nW=^Z;(s)u(s-)(i,.;i,.-|l,5,), (8) 

SS 

with its explicit form 

where the superscript (1) denotes the spin triplet state and ep is the polarization vector of the spin 1 meson. The 
Lorentz-invariant Mandelstam variables are defined by 

S = (pi +P2f = {P3+P4f, 

t = (pi -Psf = (P2 -Pif, 

U = {Pl-Pif = {J>2-P3)\ (10) 
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and they satisfy 



s + t 



P2 = 4R2 ^ 4(^2 ^ 



(11) 



Furthermore, the covariant spinors are normaUzed relativistically as uu — —vv = 2mc- 

Let Ai be short for the amphtude M.{g + g — > J/tp + g) in Eq. ([7|), and it can be expanded in powers of v or 
equivalently \q\. That is 



M 



1 8"^ MP 
e,M^ = ^P(-M''(0) + 2'^"/^ 



g=0 



) + 0(/), 



(12) 



where high order terms in four momentum q have been omitted. Terms of odd powers in q vanish because the heavy 
quark pair is in an S wave configuration. Note that the polarization vector ep also depends on q, but it only has 
even powers of four momentum q, and their expressions may be found e.g. in the Appendix of Ref. [35| . Therefore 
expansion on q^ of ep can be carried out after amplitude squaring. The following substitute is adopted: 



(13) 



This substitute should be understood to hold in the integration over relative momentum q and in the S wave case. 
Here, |q|^ can be identified as[33l[36j 



|(0|xt<T(-iD)^V>|V>(3gW))| _ (0|pJ/^(34i])|0) 
|(0|xt<T^|^(35™))| (010^/^^(35^)10) 
Then the amplitude squared defined in Eq. ([6]) up to 0{v'^) are: 



[l + Oiv% 



(14) 



J2\m\'^m'>{o)m'*{o)J2^p4 



(A 



Q/3 



r)2 A//>* 



g=0 



E 



(15) 

The heavy quark and antiquark are taken to be on shell, which means that P ■ q = 0, and then gauge invariance is 
maintained. The polarization sum in Eq. (|15p is 



-9px 



PpPx 

p2 



(16) 



It is clearly seen that the polarization sum above only contains even order powers of four momentum q, therefore 
it will make contribution to relativistic correction at 0{v'^) in the first term on the right hand side of Eq. (jlSp 
when contraction over indices p and A is carried out. However, since the second term on the right hand side of 
Eq. (|15p already has term proportional to q^, i.e. the four momentum q can be set to zero throughout the index 
contraction. Then we have 



(17) 



where A and B are independent of By comparing Eq.(|3]) and Eq.®, we obtain the short distance coefficients 
shown explicitly below. The leading order one is 



1 



1111 



167rs2 64 4 2iVc 3 
1 1111 



A 



{A-KasfbUOmclm (s^ 
2n2 



167rs2 64 4 2iVc 3 

+ (s^ +ts + t^y]/[9 {s - Amiy {t ~ 4miy {s + tf], 
and the relativistic correction term is 
G(35W) 



ts + ml -A {2s^ + 'its^ + 2,t^s 

2 



3^ ml 



(18) 



1 



ml 



167rs2 64 4 2Nc 3 
1 1111 



167rs2 64 4 2Nc 3 



(47ra,)3(-2560)[2048 (35^ + 2ts + it^) m 
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256 (55^ - 2ts^ - 2t^s + 5^^) ml 



-320 (35 



lOis^ 



lOt^s^ 



lOi^s + 3^"*) + 16 (21s^ 



63ts'* + SBt^s^ 



mt^s + 21i5) mt 



-4 Ts" + ms"" + 23t^s 



2 „4 



28t3s3 + 23t''s2 + m^s + 7t^) ml - st{s + t) (s^ +ts + t^f]/ 



[27toc (4m^ - s)-^ (4m^ - t)'' (s + t)3 



(19) 



5 



Each of the factors has its own origin: l/167rs^ is proportional to the inverse square of the M0ller's invariant flux 
factor, 1/64 and 1/4 are the color average and spin average of initial two gluons respectively, l/2Nc comes from the 
color-singlet long distance matrix element definition in Eq. ^ with Nc—i^ 1/3 is the spin average for total spin 
J = 1 states and (47ras)'^ quantifies the coupling in the QCD interaction vertices. Furthermore the variable u has 
been expressed in terms of s and t through Eq. (fTT|) . To verify our results, we find that those in Ref. |29|] discussed for 
J /ij: photoproduction are consistent with ours under replacement (47ra)e^ — > (47ras), and the result in Ref.js^] agrees 
with ours at leading order after performing the polarization summation. 



B. Nonperturbative Long Distance Matrix Elements 

The long distance matrix elements may be determined by potential model [1^ [11] or lattice calculations [111 , and 
by phenomenological extraction from experimental data. Here we first extract the decay matrix elements from 
experimental data. Up to NLO QCD and relativistic corrections, decay widths of the color-singlet J/%p to light 
hadrons (LH) and e+e~ can be expressed analytically as follows [33j: 

mi ml 
where the short distance coefficients are 1331 



^..f^W) = ~ ' , '^ ali2m.)[l + (-9.46C. + 4.13C. - 1.1617V,)^] 

iV^ lO TT 

Nf 

6 TT 

Ge+e-f^?') = (21) 

Then, the production matrix elements can be related to the decay matrix elements through vacuum saturation 
approximation 

(0|O^/'^(34'')|0) = {2J +l){H\0^'^[^S^^^)\H)[l + 0{v^)] = i{H\0-^l'^{^S^^^)\H)[l + 0[v% (22) 

Combining the above equations and the experimental data in [11], i.e., T[J/il} LH] — 81.7 KeV and V[J/%1> 
e^er\ = 5.55 KeV and excluding the NLO QCD radiative corrections in pT|) . we get the solutions accurate at leading 
order in 

^0|c)J/V'(35.[i])|o) = o.868GeV3, (0|7"^/'^(3S'P')|0) = 0.190GeV^ (23) 

And the enhanced matrix elements accurate up to NLO in can be obtained by including NLO QCD radiative 
corrections in (PT|) 

^O|0,//^(35.[i])|Q^ = lMGeV\ (0|p-'/'^(3S'™)|0) = 0.32006^^ (24) 

The strong coupling constant evaluated at the charm quark mass scale is as{2mc) — 0.250 for mc = 1.5GcV. The 
other input parameters are chosen as follows: the QCD scale parameter Aqcd = 392MeV, the number of quarks with 
mass less than the energy scale mc is Nf — 3, color factor Cp — 4/3 and Ca — 3, the electric charge of charm quark 
is eg ~ 2/3, Qi are the electric charges of the light quarks and fine structure constant a,,, = 1/137. Our numerical 

values for the production matrix elements (0|O'^/'^(3S'[^')|0) and (0|7''^/'^(^S'P')|0) are accurate up to NLO in V with 
uncertainties due to experimental errors and higher order corrections. 
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C. Cross Sections for p + p{p) J/tp + X and Phase Space Integration 

Based on the results obtained for the subprocess g + g J/tp + g we further calculate the production cross sections 
and relativistic corrections in the process p + p{p) J/V' + X , which involves hadrons as the initial states. In order 
to get the cross sections at the hadron level, the partonic cross section defined in Eq.® has to be convoluted with 
the parton distribution function (PDF) fg/p{x±), where x± denotes the fraction of the proton or antiproton beam 
energy carried by the gluons. 

We will work in the pp center-of-mass (CM) frame and denote the pp energy by \/S, the rapidity of J/iIj by yc, 
and that of the gluon jet by yo- The differential cross section of p + p{p) — > J/ip + X can be written as|39l| 

d^a{p+p{p) ^ J/^ + X) r , . r , , da{g + g^ J/i, + g) 

^p2^y^^y^ - x+/,/,(x+)x_/,/,(,)(x_) , (25) 

where 

_ m^exp{±yc) + 'm^exp{±yD) 



with the transverse mass mf'^ ~ ym^jy+p^, the meson mass nic — 2mc, the gluon mass niD — and the 
transverse momentum pt. The Mandelstam variables can be expressed in terms of pt, yc and yo 

s — X-^x — iS*, 

t = -pj-mfmj^exp{yD-yc), 

u = -p1 - m^m'^expiyc - yo)- (27) 

The accessible phase space puts kinetic constraints on variables pt, yc and yu for a fixed value of two colliding hadron 
center-of-mass energy 



< „ < y KS.^h.'^'n) 

S + — mj~i 



\yc\ < Arcosh 



c 



, - m^exp{-yc) . . , y/S - mf expjyc) 
-In g < < In 5 , (28) 

where A(a;, y, z) = x'^ + y'^ + z'^ — 2{xy + yz + zx). The distribution over pt of differential cross section can be obtained 
after phase space integration. 



III. NUMERICAL RESULTS AND ANALYSIS 

The CTEQ6 PDFs[40l| are used in our numerical calculation. We present the distribution of J/ip production 
differential cross section da/dpt over pt at the Tevatron with ^/S = 1.96 TeV and at the LHC with = 14 TeV in 
FigUl Figl3]and FigHl FiglS) The solid line represents the distribution at leading order in 0{v'^), and the dotted line 
describes the relativistic correction at next-to-leading order in 0{v'^) (excluding the leading order result). The long 
distance matrix elements are accurate up to leading order in as from Eq. (|23p or next-to-leading order in as from Eq. 
([M)) . The variable pt is set to be from 5 GeV to 30 GeV (50 GeV) for the Tevatron (LHC), and the distributions are 
depicted in logarithm unit along the vertical axis. All curves decrease rather rapidly as the transverse momentum pt 
increases, and the leading order da/dpt behavior is not changed by the relativistic corrections. It can be seen that the 
ratio of relativistic correction to leading order term is 1% or so, and less than 2%, which is insignificant and negligible. 

The tiny effect of relativistic corrections is partly due to the smallness of the short distance coefficient correction. 
In fact, the ratio of NLO short distance coefficient to LO one from Eq. ((19)) and Eq. (fT8|) can be expanded as a series 
of the small quantity mj,, as compared with ^/s, and this series reduces to a fixed small number ^ if only the leading 
order term is kept, i.e., 

^Mll^l as ^^O^^O (29) 



7 




1E-6 r 



1E-7 r 



1E-8 I . 1 . 1 . 1 . 1 . 1 

5 10 15 20 25 30 

P. (GeV) 

FIG. 2: The pt distribution of da{p + p J/'tp + X)/dpt at the Tevatron with %/S = 1.96 TeV. The 0{v°) and 0(v^) resuhs 
are represented by the sohd and dotted hnes respectively. 
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FIG. 3: The pt distribution of da{p+p J/ip + X)/dpt (with enhanced matrix elements) at the Tevatron with y/S = 1.96 TeV. 
The 0{v^^) and 0{v'^) results are represented by the solid and dotted hnes respectively. 



Together with the suppression from long distance matrix elements, the tiny effect of relativistic corrections can be 
accounted for. Our results for relativistic corrections in the process p + p{p) — > J/ip + X are similar to that in the J/tp 
photoproduction process discussed in Ref.f^o']. These results may indicate that the nonrelativistic approximation in 
NRQCD is good for charmonium production at high energy collisions, and relativistic corrections arc not important. 
This is in contrast to the case of double charmonium production in e+e^ annihilation at B factories, where relativistic 
corrections may be significant. 

IV. SUMMARY 

In this paper, relativistic corrections to the color-singlet J/ijj hadroproduction at the Tevatron and LHC are calcu- 
lated up to 0{v'^) in the framework of NRQCD factorization approach. The perturbative short distance coefficients 
are obtained by matching full QCD differential cross section with NRQCD effective field theory calculation for the 
subprocess g + g J/ip + g- The nonperturbative long distance matrix elements are extracted from experimental 
data for J/i/j hadronic and leptonic decay widths up to 0{v'^) with an approximate relation between the production 
matrix elements and decay matrix elements. Using the CTEQ6 parton distribution functions, we then calculate the 
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FIG. 4: The pt distribution of da{p + J/ip + X)/dpt at the LHC with ^ = 14 TeV. The 0{v°) and Oiy'^) results are 
represented by the solid and dotted lines respectively. 



> 



T3 



■LO(enhanced) 
Correction(enhanced) 
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FIG. 5: The pt distribution of d(T(p + p ^ J/-i/> + X)/dpt (with enhanced matrix elements) at the LHC with \/S — 14 TeV. 
The 0{v'^) and 0{v'^) results are represented by the solid and dotted lines respectively. 



LO production cross sections and relativistic corrections for the process p+p{p) — > J/ip + X at the Tevatron and LHC. 
We find that the 0{v'^) relativistic corrections to the differential cross sections over a wide range of large transverse 
momentum pt are tiny and negligible, only at a level of about 1%. The tiny effect of relativistic corrections is due 
to the smallness of the short distance coefficient correction and the suppression from long distance matrix elements. 
These results may indicate that the nonrelativistic approximation in NRQCD is good for charmonium production at 
high energy hadron-hadron collisions, and relativistic corrections can not offer much help to resolve the large discrep- 
ancy between leading order prediction and experimental data for J/'ijj production at the Tevatron. Other mechanisms 
such as those suggested in jJH, \^ may need to be considered, aside from higher order QCD contributions. 
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